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This paper presents a unified semiclassical framework for five regimes of quantum fidelity decay 
and conjectures a new universal regime. The theory is based solely on the statistics of actions in 
the dephasing representation. Counterintuitively, in this representation, all of the decay is due to 
interference and none due to the decay of classical overlaps. Both rigorous and numerical support 
of the theory is provided. 
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While two classical trajectories with sUghtly different 
initial conditions will exponentially diverge in a generic 
system, overlap of two quantum states evolved with the 
same Hamiltonian H remains constant for all times. The 
situation changes if we consider the sensitivity to pertur- 
bations of the Hamiltonian. In a generic system, two tra- 
jectories with the same initial conditions, but propagated 
by slightly different Hamiltonians i?" and H'^ will also ex- 
ponentially diverge. The quantum-mechanical version of 
this question was posed only very recently by Peres P| 
and has been extensively studied in the past three years 

Jk ^, ^ ^ i A ^ Jii, ii, 12* ill Q El E E3 El El 

1201 especially due to its relevance to quantum computa- 
tion '2l'|. The quantum- mechanical sensitivity to pertur- 
bation is measured in terms of quantum fidelity M (t) : 
the overlap at time t of two initially identical states \ip), 
that were evolved by slightly different Hamiltonians i?" 
and iJ' — -I- eV (e controls perturbation strength). 
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Recent investigations of the temporal decay of fidelity 
discovered a plethora of regimes in integrable and chaotic 
systems. The decay can be Gaussian ^ |1|, exponen- 
tial 0, H 01 or superexponential El in chaotic systems, 
and Gaussian 01 or algebraic H in integrable systems. 
There exist actually two qualitatively different exponen- 
tial regimes in the chaotic systems: the Fermi-Golden- 
Rule regime where the decay rate is proportional to 
0, II, and the Lyapunov regime with a perturbation- 
independent decay rate 0. Many authors are unaware 
of the difference between the perturbative Gaussian decay 
in chaotic systems (which occurs only after the Heisen- 
berg time |6|) and the Gaussian decay in integrable sys- 
tems (which occurs much before the Heisenberg time 
Generic dynamical systems are mixed: their phase space 
consists of both invariant tori and chaotic regions, and so 
in general fidelity decays by a mixture of the five regimes. 

Various regimes have been qualitatively explained 
by different quantum, statistical, or semiclassical tech- 
niques, but no single technique described more than two 
regimes. Moreover, semiclassical explanations always 
made an assumption about the decay of classical over- 



laps. This paper presents a unified semiclassical frame- 
work [Eq. 0] for five known regimes, using exclusively 
the statistics of actions in the dephasing representation, 
rendering extra assumptions unnecessary. 

In Ref.|3|, a simple uniform (i.e., free of singularities) 
semiclassical (SC) expression was derived for fidelity of 
initial position states and Gaussian wave packets. In 
Ref.'s), this expression was rigorously justified by the 
shadowing theorem, and generalized to any Wigner dis- 
tributions, describing both pure and mixed states. Due 
to the lack of a SC dynamical prefactor, the approxima- 
tion was called the dephasing representation (DR), 
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Here x = (r, p) denotes the phase space coordinates, r' 
and p' are the initial position and momentum, and 
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is the difference at time t of actions of the perturbed and 
unperturbed Hamiltonian along the unperturbed trajec- 
tory which is equal to the negative of the integral of the 
perturbation along the unperturbed trajectory. 

Although Eq. works for any Wigner distribution, 
let us consider first the most random state pw = 
fidelity of which best approximates fidelity averaged over 
initial states. Here O stands for the phase space volume. 
The central result of this paper is that the five regimes 
occurring before the Heisenberg time can be described by 
a single formula. 



Mdr (t) 



j2d 



a;_ exp 



1 

'2tf 



where AS" = AS'(x',t), AS"' = AS'(x",t) , x_ 
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x' — x", and the average is over x_| 



(x' + x")/2. 



All one must do is find the 
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variance ((AS'— AS") ). The two criteria that yield four 
different regimes are first, whether AS' and AS"' are cor- 
related, and second, whether the dynamics is chaotic or 
quasi-integrable. Table shows the variance in the four 
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TABLE I: Variance ((AS" - A5"')^). 



Dynamics 
Chaotic Quasi-integrable 



A5" and AS": 
Uncorrelated 
Correlated 
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TABLE II: Fidelity M{t). fip are volumes of the most 
unstable direction and of the momentum space; (3, 7 are in- 
dependent of e and t. 
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different cases. If one substitutes an appropriate expres- 
sion for variance from Table into Eq. ^ and performs 
the trivial integral, one obtains Table UTI of temporal de- 
cay of fidelity in the four regimes. Table ITTll associates 
these decays with the terminology used in literature. 

To obtain the universal expression Q , one starts from 
the overlap Q squared for random states pw = , 



Mdr [t) 



1 



S''x' I d'^'x" exp 



(A5'-AS"' 



(4) 



Changing variables to xj- and averaging over x+ gives 
1 



Mdr, it) = 



d?'^x_ ( exp 



(AS"-AS"') 



If we could isolate each regime, the difference AS"— AS"' 
would be generally Gaussian distributed (even if AS" itself 
were not). Since (AS"-AS") = 0, when AS"-AS"' is 
Gaussian distributed, we get the universal Eq. (Q. 

What remains is deriving entries in Table AS' 
and AS" are uncorrelated if the initial coordinates are 
different enough [x^ large) or for long enough time t 
when the initial correlation is forgotten. Then the vari- 
ance of the difference is just twice the variance of each 



TABLE III: Regimes of fidelity decay. 





Dynamics 




Chaotic 


Quasi-integrable 


AS' and AS": 






Uncorrelated 


Fermi-Golden- Rule 


Gaussian 


Corr. (large t) 


Lyapunov 


Algebraic 


Corr. (small t) 


Superexponential 


(Cubic-exponential) 



term, ((AS'-AS'") ) = 2cr|g, where AS' is given by 
Eq. (0). Crucial quantity is the potential correlator 
Cv [t) ~ {V [r [t)] V [r (0)])q. In chaotic systems, or any 
systems in which this correlator asymptotically decays 
faster than 1/t, AS follows a random walk, and 
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with K = Jp°° dt Cv (t) 0, 01 ■ In quasi-integrable sys- 
tems, or systems where Cy asymptotically oscillates 
about a finite value, AS follows a ballistic motion, and 
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where Cy = Umt^tx. /g drCy (t). Analogous result 
was obtained quantum-mechanically by Prosen . 

AS" and AS" are correlated if the initial coordi- 
nates are close enough (x- small) or for short enough 
time before the initial correlation is forgotten. Since 
((AS'— AS'")^) 7^ 2(T^g, we cannot use a simplifica- 
tion as above. Nevertheless, AS"— AS" itself follows 
a generalized random walk with a time-dependent step 
(xSr (i) = r'(t) -r"(t). 
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where r(t) = (l/2)[r'(t) -I- r"(i)]. If the force-force cor- 
relator Cpit) = (Vy [r(t)] • Vy[r(0)])f2 decays faster 
than 1/t then ((AS'-AS")^) « 2De^ dr Sr (rf with 
D = dtCp (t). In chaotic systems, Sr (t) w axue"^* 
where Xu is the projection of x- onto the imstable direc- 
tion and a depends weakly on t. Then 



((A^'-A^")') 



xl/X. 



(5) In quasi-integrable systems, Sr{t) ^ p_t/m, and 



{{AS'~AS"y) w 2DeVt73m 
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Although it is safer to assume that Cpit) decays faster 
than 1/t, this assumption seems over- restrictive because 
D can be finite even if Cpit) decays as 1/t or slower, as 
long as it oscillates (e.g., consider Cpit) ^ sin t/t). 

For long enough times, the size of phase space is ir- 
relevant and the limits of the integral in Eq. Q can be 
replaced by infinity. When the results I©, Q), lO, ifTUI) 
are substituted into the universal Eq. (PJ, we obtain the 
first two rows in Table ITTl and the corresponding regimes 
in Table Jjllll The dependence on t and e agrees with 
Refs. HlllEll, 0, H, and 0, respectively. 

For short times, when even trajectories with most dis- 
tant initial conditions are still correlated, size of phase 
space comes into play, and substituting variances lO and 
(fTn|l into Eq. © gives the last row of Table lul and Ta- 
ble HHJ The first result agrees with Ref. |l5l |. the latter 
predicts a short-time cubic-exponential decay in quasi- 
integrable systems, yet to be observed. 
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FIG. 1: Variance of AS as a function of time in a a) chaotic 
and b) quasi-integrable system. Dependence is a) linear (slope 
^ 0.992) and b) quadratic (slope ?a 1.9995). 

While the DR ^ has already been tested for position 
eigenstates in chaotic 0, 0, 0| and mixed systems, 
detailed numerical tests of the more general expression 
for both pure and mixed states, in chaotic, mixed, 
and quasi-integrable systems, in the first five regimes 
from Table IIIII as well as intermediate regimes will be 
presented elsewhere Here the numerics is focused 

on the verification of Table ^ which is the only input 
required by Eq. In particular, the dependence of 

on p_ and t is checked. Many authors 
believe that it is necessary to calculate fidelity averaged 
over initial states to obtain regimes in the second row of 
Table Unl To disprove that belief let us use pure position 
states pw — ^p^S{r — R). Derivations in Eqs. Q-© 
will stay the same if we replace integrals D,^-^ J d^'^x ■ ■ ■ 

For numerics, a perturbed standard map [Sj, 

Pn+i ^ Pn + k sin q„ + e sin 2q„ (mod 2tt), 
Qn+i = qn +Pn+i (mod27r). 

was used. Here g„, Pn are the position and momen- 
tum at discrete times n, and A; is a parameter control- 
ling the transition from integrability to chaos. In the 
calculations described below, fc = 20 was used as a rep- 
resentative chaotic system, and k = 0.3 as a represen- 
tative quasi-integrable system. The initial state is a po- 
sition state Q — 0.8 tt, avoiding any problems due to 
symmetry. Between 500 and 1000 trajectories were used 
to check the statistics of actions. Effective Planck con- 
stant is h ~ l/27m where n is the size of Hilbert space. 
Specifically, n — 1000 and e = 0.003 were used for the 
chaotic example, and n ~ 100 and e = 0.005 for the 
quasi-integrable case. 

Dependence of aj^g on time, described by Eqs. ® 
and O , is verified in Fig. While part a) shows that 
in chaotic systems, a'^g grows linearly with time, part 
b) shows that in quasi-integrable systems, a'^g grows 
quadratically with time. Fig. |21 verifies the dependence 
of ((A5" — AS"')^) at a fixed time on the difference p- 
of initial momenta. Parts a) and b) show that in both 
chaotic and quasi-integrable systems, this dependence is 
quadratic for small p_ (as in Eqs. © and (|10|l ) and in- 
dependent of p_ for large p_ (as in Eqs. © and Q). 
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FIG. 2: Variance of AS' — AS" at a fixed time t = 7 as a 
function of j)_ = p' —p" in a a) chaotic and b) quasi-integrable 
system. In both cases the dependence is first quadratic and 
then independent of p' —p". Fitted slopes: a) 1.998 and 0.01, 
b) 2.000 and 0.001. 
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FIG. 3: Variance of AS' — AS" as a function of time for fixed 
p' —p" — 10~^ in a a) chaotic and b) quasi-integrable system. 
The dependence is: a) first exponential, then linear, and b) 
first cubic (slope « 3.046), then quadratic (slope ~ 1.994). In 
b) , dots represent the numerical variance averaged over a time 
interval {t/2,t) since the variance itself has large oscillations 
about the averaged value. 

The transition occurs for p_ such that by the time i, two 
trajectories with initial distance p_ completely lose their 
correlation. Dependence of ((AS" — AS"')^) on time for a 
given difference of initial momenta is confirmed in Fig. 
Part a) shows that in chaotic systems this dependence 
is first exponential, as in Eq. ©, and later linear, as in 
Eq. © (hard to see here, but can be seen in a log-log plot 
such as in Fig. © . Part b) shows that in quasi-integrable 
systems, this dependence is first cubic, as in Eq. H10|l . 
and then quadratic, as in Eq. ©. The transition occurs 
at time t when two trajectories with initial distance p„ 
completely lose their correlation. 

Intuitively, one would think that the decay of quan- 
tum fidelity would have two components: the decay of 
classical overlaps (classical fidelity) and the decay due to 
dephasing (destructive interference). This approach was 
taken in all previous fidelity literature. Contrary to that, 
in the DR, all of the decay is due to interference. This 
can be best seen in the algebraic decay: while in Ref. |^ , 
the decay of classical overlaps, oc t~'^, and the decay due 
to dephasing, oc t"''/^, together give the overall fidelity 
decay, oc t~^'^/'^ , in the present approach, the same overall 
decay t~3d/2 g^tirely due to dephasing in DR. 

I would like to point out the rigorousness of going from 
((AS"-AS"')^) to M (t) in the DR approach: here the in- 
tegration is done analytically and correctly. Similar ap- 
proach based on the statistics of action differences (how- 
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ever, in the final position representation where the inter- 
ference accounts only for a part of the decay) was used 
in literature to derive the Lryapunov and algebraic 
decay. Both Refs. |3| and j^j provide a long derivation in 
which there appears an integral of the form 



J d'^rJ2 |dct(ap^./9r) 



in particular, there is a second power of a certain Jaco- 
bian (Van Vleck determinant). In both Refs. 4, 9], a 
suspicious change of variables is employed, in which one 
power of the Jacobian is used to change variables from r 
to p' and the other power of the Jacobian is replaced by 
its estimate to get 



d'^p' |det(ap7(9r) 



estimate 



While both Refs. obtain the right overall behavior at the 
end (because they guess the behavior of the Jacobian), 
this step makes the lengthy mathematical derivation non- 
rigorous. 

Originally, dephasing representation JQ) was used only 
in chaotic systems 2^. Besides generalizing DR to arbi- 
trary Wigner functions, the main result of Ref. .3| was 
showing that DR is valid for finite times also in mixed 
and quasi-integrable systems. The mathematical cor- 
nerstone of the approximation is the shadowing theo- 
rem [23, lU I23. 125L Eif which guarantees the existence 
of an unperturbed trajectory with slightly different ini- 
tial conditions that is near (shadows) a perturbed tra- 
jectory up to time ts- In uniformly hyperbolic systems, 
time ts is infinite, but for a much larger class of sys- 
tems, shadowing works for finite times tg. In its origi- 
nal form [23, 123, m, 113, shadowing theorem justi- 
fies computer-generated trajectories in chaotic systems 
where perturbation is a random noise due to the round- 
off errors. The relation with deterministic Hamiltonian 
perturbations was rigorously shown in Ref.Q- 

Dephasing representation ^ appears too simple to be 
true because Eq. (0) has the form of first order pertur- 
bation approximation and appears to consider only the 
interference of the "diagonal" terms. The DR seems to 
require that the change of the trajectory due to pertur- 
bation be very small. While some authors believe that, 
this is not true at all. Other authors believe that the ap- 
parent first-order perturbation approximation is a conve- 
nience and that it could be improved by considering the 
exact trajectories. That is not true either. The reason 
why DR works is the shadowing theorem as explained in 
Ref.Q- While a perturbed trajectory with the same ini- 
tial condition can change enormously, there exists a per- 
turbed trajectory with slightly different initial condition 
that changes very little. Trying to improve the approxi- 
mation by relieving either apparent ingredient ( "diagonal 
approximation" or "first order perturbation approxima- 



tion") not only does not improve the numerical results, 
but actually gives much worse results. 

To conclude, a simple unified framework was presented 
for the temporal decay of quantum fidelity. It is based 
solely on the statistics of actions in the DR, through 
Eq. and can describe five known and one new univer- 
sal regime. The DR ([Ql 0, Q is even more general and 
can represent quantum fidelity in a mixture of regimes as 
well as in other non-universal regimes 17]. 

The author wishes to thank the Mathematical Sciences 
Research Institute for financial support. 
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